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Abstract 

We present exact results for the density profile of the one dimen- 
sional array of hard spheres of diameter D and mass m under gravity 
g. For a strictly one dimensional system, the liquid-solid transition 
occurs at zero temperature, because the close- packed density, (pc, is 
one. However, if we relax this condition slightly such that (j)c = 1 — S, 
we find a series of critical temperatures, Ti*^ = ■mgD{N + 1 — i) / fig 
with iio = 1/6 — 1, at which the z-th particle undergoes the liquid- 
solid transition. The functional form of the onset temperature, tJ^^ = 
mgDN/ fio, is consistent with the previous result [Physica A 271, 192 
(1999)] obtained by the Enskog equation. We also show that the in- 
crease in the center of mass is linear in T before the transition, but it 
becomes quadratic in T after the transition because of the formation 
of solid near the bottom. 

PACS: 64.70Dv; 51.30+i; 45.70vn 



I. Introduction 

In the previous paper [1], the author proposed that the hard sphere gas 
undergoes the condensation transition under gravity g, and identified the 
transition temperature, Tc, as the point at which the Enskog equation [2] fails 
to conserve the total number of particles. Based on the fact that hard spheres 
cannot be compressed beyond the close-packed density, it was suggested [1] 
and confirmed [3] that the missing particles should condense from the bottom 
and form a solid below T^, and its fraction in a solid regime at a temperature 



1 



T < Tc was predicted to be 1 — T/Tc. The transition temperature Tc was 
determined as 

Te = mgDii/no (1) 

where m, D are the mass and the diameter of a hard sphere, /i is the layer 
number of the system, and Ho is a constant, which depends on the level of 
approximations in truncating the BBGKY heirarchy [4], or perhaps in em- 
ploying the density functional theory [5]. Hence, the value obtained by the 
Enskog theory [1], or more precisely by the Enskog pressure, may be close 
to the real value, but not precise. For example, if one uses the pressure form 
suggested by Percus-Yevick [6], this constant will be slightly different. For a 
one dimensional lattice gas [7,8], it can be shown that fio = —ln{a/{l — a)) 
with a = exp{ — 14). The crucial point, however, is that the scaling form 
of the transition temperature (Eq.(l)) should survive in all approximations. 
The purpose of this paper is to demonstrate this point by exactly solving 
the one dimensional hard sphere problem under gravity. For a strictly one 
dimensional system, the condensation transition occurs at zero temperature, 
because the close-packed density, (f)^, is one. However, some useful informa- 
tion may be extracted from Id results if we relax this condition slightly such 
that (pc — i — S. Then, the transtion occurs at a finite temperature. Even 
though the fluctuations in 3d are very different from Id system, we will show 
in this paper that results obtained in this way appear to be relevant to the 
real physical system. Perhaps, one may view such a Id system as a coarse 
grained mean field system of the real three dimensional hard sphere system. 
We will obtain the exact transition temperature, T^, and check its functional 
form against Eq.(l). We also determine the series of transition temperatures, 
T^^\ at which the i-th layer undergoes the condensation transition. We fur- 
ther show how a sharp departure in the center of mass statistics shows up 
before and after the transition. Before the transition, the increase in the 
center of mass is linear in T, while after the transtion it is quadratic in T, 
because of the formation of solid near the bottom, which is a characteristic 
of Fermi systems [9,10]. 

II. Condensation of One Dimensional Hard Sphere Gas Under 
Gravity 

Consider a collection of hard spheres of finite radius R (or diameter D = 
2R) in a one dimensional tube with the top open. Let the mass of the i-th 
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particle be rrii. We assume that each hard sphere is in thermal equilibrium 
with a heat reservoir at a temperature T. The system we have in mind is 
the one used in usual Molecular Dynamics simulations, where each particle 
is kicked periodically by Gaussian noise so that the average kinetic energy 
of each particle, m < v"^ > /2 = T. We ignore the pressure due to the 
reservoir. In such a case, since the kinetics is separated out, we only consider 
the configurational integral in computing the partition function, Z^r, of the 
N particle assembly: 

Zn = / dzi / / dzNexp{—P'g{miZi + + rnNZN)) (2) 

JR Jzi+2R Jzm-i+2R 

with P' — 1/T. The hard sphere gas without gravity has been studied and 
is known as the Tonk gas [11]. The integral in (2) involves the exponential 
functions and thus can be carried out exactly to yield, 

]^ g-2/3'pmiv-R g-2/9'g(mjv+mjv_i)R g-2/3'p(mi+m2+...+mjv)ii 

{(5'gY rriN {rriN + m^v-i) (mi + m2 + + m^v) ^"^^ 

We now compute average quantities. First, the average position of the i-th 
particle < > is given by: 

< Zi >= = 2z - l)R + -z, 4 

(5'g drrii g 

where 

If all the masses are the same, i.e. = m, then this reduces to 

<Zi> /D = {z- 1/2) + ^— ^ : (6) 

Note that the first term, Zi{0)/D — i — 1/2, results from the close-packing 
at the ground state T = and the second term represents the thermal ex- 
pansion. Note also that Y.i < Zi > / D = + TN/mgD. The dimen- 
sionless thermal expansion defined as, Zi — (< ^Zi > / D){mgD /T) with 
< /S.Zi >—< Zi{t) > —Zi{0), is independent of the temperature. For example, 

z^ = l/N, Z2 = l/N+1/{N-1), ZN = l/N+l/{N-l) + .... + 1/2 + 1 
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The dimensionless mean expansion per particle is precisely given by the 
thermal energy injected into the system: 



i i 1=1 

(7a) 

The change in the center of mass due to the thermal expansion is linear 
in T: 

1 ^ 

< z{T) >= ^ E(< ^^iT) > -ZiiO)) T/mg (76) 
i=i 

Wc now compute the density profile, p{zi), as a function of position, Zi. 
Define the dimensionless density, (f){zi) = p{zi)/pc with = 1/D. Then, 
since p{zi)Azi = Ai, we find (p{zi) = -^/Pc and its discrete version becomes: 

*W = Vli + i^^] (8) 

where we have used the relation: 

rx=N—i-\-l 

l/N+l/{N-l) + ... + l/{N-i + l) ^- dx/x^ -ln{l-{i-l)/N) 

Jx=N 

(9) 

and we have redefined the dimensionless temperature (5 — mgD/T and the 
dimensionless coordinate, yj =< Zi> / D: 

y^-i'-m + \tjrz^, (10) 

Note that / p{zi)dzi = J ((){yi)dyi = N(f)c with 0c the close-packed density. 
The density at the bottom layer, 0O, is given by Eq.(8) with i = 1, i.e., 
4>o= + V/^^]- -^o^ ^ strictly one dimensional system, the closed packed 
density, 0c = 1, and thus by setting, 0o = 0c = 1, we find that the crystal- 
lization occurs at zero temperature in one dimension. 

In order to extract some useful information from one dimensional results, 
and make them relevant to higher dimension, we assume that the close- 
packed density is slightly below one by a small amount, 0<57^0<<1, 
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i.e., 0c = I — S. What we have in mind is a coarse grained three dimensional 
system, for which each column may interact weakly. In fact, we have found 
that such a system can be realized in the Molecular Dynamics simulations 
if the system is initially arranged in a two dimensional square lattice with a 
little space between the columns. In such particles in each column 

do not mix, and the square structure is maintained [18]. Such a model could 
be understood as a coarse grained mean field model in the spirit of rcf.[7]. 
Certainly, the fluctuations in 3d are very different from those in Id, and 
thus it may be objectionable to extent the results of Id to 3d. Nevertheless, 
the results of the Id obtained this way with regard to the existence of the 
condensation temperature, and perhaps the existence of the discrete jump in 
the condensation process may survive in high d, as will be shown shortly. 

Now, if we let 0c = 1 — then one can easily find the onset of the 
condensation temperature, Tc, at which the first layer becomes crystallized: 



where /i — N is the initial layer number (or the Fermi energy [10]), and the 
constant is given by 



Note that Eq.(ll) has the same funtional form as (1). One may relate S 
to the critical pressure, PcD^ — /loTc/ D at which the crystalization occurs. 
From the force balance equation, we find the pressure at the bottom: 



The factor D'^ was introduced to effectively model the three dimensional 
system. By equating P(0) to the critical pressure, Pc, we again find the 
transition temperature, = mgDfi/ Ref. [7] identifies the critical pressure 
as: Pc = UTJD^. Hence, we find: /x„ = 1/5 - 1 = 14.937, and 5 !^ 
0.06275. After the first layer becomes cystallized at Tc, the density profile 
above the second layer is given by Eq.(8) with N replaced hy N — 1, and 
i = 1, ...,N — 1. This is effectively equivalent to shifting the origin from the 
first to the second layer. The second layer, which has now become the origin, 
becomes crystallized at the second critical temperature, T^^^: (po{T^'^^) = 1—6. 
The process continues, and we can find a series of critical temperatures, T^^\ 



Tc = mgDn/iio 



(11) 




(12) 




(13) 
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at which the i-th layer in the original labeling becomes crystallized: 



^ mgD{N+l-t) 



So, all the particles are crystalized at T = Tj:^^ = mgD/no, which is not the 
absolute zero, because 5 7^ 0. Note also that the crystallization of each layer 
proceeds with a discreet temperature jump, AT = T^(*+^) — T^(*) = mgD/ jio- 
Hence, the heat release, or the latent heat, Q, resulting from the formation 
of one solid layer is Q = AT = mgD/no- Biben et al [12] investigated 
the density profile of a hard sphere suspension in a gravitational field using 
Monte Carlo simulations, and reported that for A = mgD/T < 2.5, the 
system is a strongly perturbed fluid, while at A fa 2.75 the first two layers 
form a crystal, and the formation of third and fourth layer crystals occurs in 
a discontinuous manner between A = 2.5 and A = 2.75. Setting mgD/Tc — 
2.5, and T^. = mgD/fio, we find fig = 2.5 and 6 ~ 0.2857, and the gap 
AT = mgD/no ~ 0.25. Such findings do not seem to be inconsistent with 
the results presented above. 

We now examine the center of mass statistics below the condensation 
point Tc = T^^\ At a given temperature tJ*"*"^) < T < TJ*\ what is the 
fraction of particles in a condensed regime? At this temperature, particles 
up to the i-th layer are condensed. Then, the fraction of particles in the 
condensed regime, = i/N, which is termed the Fermi surface [1], is given 
by a simple manipulation of identities: 

C^ = i/Ar=l-[^] = 1-T/Te (15) 

where we used, T^'+^^T^ = [mgD{N - i)/ n^]/mgDN/ n^] = {N -i)/N = 
TIT,. 

Now, the dimensionless center of mass, < y(T) >=< 2;(T) > /D,\s given 

by: 

/•CO roo 

< y{T) >= / dyy<P{y)/ / dy<P{y) = h/h (16) 

Jo Jo 

where the integral now splits into two due to the formation of a solid below 
(f- More precisely, 

h = / (l>cdy + (t>{y- CF)dy = 0cCf + {N - Cf)0c = Ncf), (17) 

Jo J(^p 
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We need some manipulation in computing the denominator I2. To this end, 
we again spht I2 into two integrals: one for the solid regime, which is es- 
sentially a rectangle, and the other for the fluid regime, where the density 
profile is given by (8) but with N replaced hy N' — N — (p- Hence, 

h = / y(t>cdy + / yckiy - CF)dy = 0cCf/2 + CF(t>c{N - Cf) + J (18) 

Jo J C,F 

where 

.00 N-^F A N-Cf 

/ dyyct>{y)= ^ y^ct>,{^),=, = ^ z,/D (19) 

But Ef Zj = D[N'^/2 + TN'/mgD]{E(i.{Q)). Hence, with N' = N - Cf = 
NT/T^, we find: 

J = [N^T^/2Tj + {N^T^ / iioT^)] = N^k(T/%f (20) 

where A = [1/2 + I///0]. Note that the increase in the center of mass is 
quadratic in T, namely: 

< A^(T) >=< z{T) > -ND/2 = aND{T/Tcf (21) 

with a = [{2 + 5{1 — 5)]/[2(l — 5)^], which is a characeristic of Fermi systems 
[10]. 

In passing, we make the following remarks. In ref. [7], an attempt was 

made to derive the condensation point for the lattice gas, which is again 
consistent with the form given by Eq.(l). While the lattice gas may capture 
some of the essence of the hard sphere systems, it is important to recognize 
that the logarithmic singularity in the pressure of the lattice gas [7,13] is far 
different from the power law singularity of the real hard sphere gas [14] . Fi- 
nally, the relevance of the present study to granular materials [15]: Granular 
materials are macroscopic particles, and the parameter A = mgD/T ^ 10^^ 
is an astronomical number, if one uses a usual temperature. Hence, the tem- 
perature T of the hard sphere gas should be interpreted differently. One way 
to relate this temperature to the vibrational strength of the granular bed is 
to compare the kinetic expansion of the granular bed to the thermal expan- 
sion of the hard spheres, as was done in ref. [10]. If we denote by h{r) the 
jump height of a single ball in the vibrating bed of the vibrational strength 
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r = AuP I g with A and a; the amphtude and the frequency of the vibration, 
then we may set 



< Az(r) >= aND{T/T^f = h{T) (22) 

from which we can find the relation between the thermal temperature T of 
the hard spheres and the vibrational strength V: 



7; ^ '\DN2 + 25{1-5) 
or, equivalently: 



" ' ' (23) 



^^^^ (24) 



mg 2 + 5(1 - 5) 

We point out that for granular materials excited by vibration in a two dimen- 
sional container, A = mgD/Tc was determined by fitting the density profile 
of rcf.[9] by the Enskog profile. The estimated value was A 4.926 [3], 
and the dimcnsioncss temperature of the vibrating bed was: T/Tc = 0.663. 
However, we point out that we have not taken into account (a) the internal 
degrees of freedom [16] of the macroscopic particles, such as rotation, and 
(b) the inelastic collisions, which may lead to an interesting clustering insta- 
bility [17]. Hence, one has to be somewhat careful in extending the results 
of elastic hard spheres to granular materials. 
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